Evidence for multi-band strongly coupled superconductivity in SmFeAsOo.sFo.2 single 
crystals by high-field vortex torque magnet ometry 
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To probe manifestations of multiband superconductivity in oxypnictides, we measured the angular 
dependence of magnetic torque t(9) in the mixed state of SmOo.8Fo.2FeAs single crystals as functions 
of temperature T and high magnetic field H up to 30 T. We show that the effective mass anisotropy 
parameter 7 extracted from t(9), can be greatly overestimated if the strong paramagnetism of Sm 
or Fe ions is not properly taken into account. The correctly extracted 7 depends on both T and 
H, saturating at 7 ~ 9 at lower temperatures. Neither the London penetration depth nor the 
superfluid density is affected by high fields fields up to the upper critical field. Our results indicate 
two strongly-coupled superconducting gaps of nearly equal magnitudes. 

PACS numbers: 74.25.-q, 74.25. Ha, 74.25.Op, 74.70.Dd 



The recently discovered superconducting oxypnictides 
[3, Q have similarities with the high T c cuprates, such as 
the emergence of superconductivity upon doping a parent 
antiferromagnetic compound @, y, 5 . Several theoreti- 
cal models [fa, @| suggest unconventional superconducting 
pairing, while the Andreev spectroscopy [7|, penetration 
depth [8J, and photoemission measurements [9|indicate 
nodless s-wavc pairing symmetry. Experiments [9l. fiol. fill] 
have found evidence for multi-gap superconductivity, in 
agreement with theoretical predictions [Bj]. 

The comparatively high T c values and extremely high 
upper critical fields H C 2 of the oxypnictides [l(J E3] in- 
dicate promising prospects for technological applications 
if, unlike the layered cuprates, a sizeable vortex liquid 
region responsible for dissipativc flux flow does not dom- 
inate their temperature-magnetic field (T — H) phase 
diagram. It is therefore important to reveal the true 
behavior of the anisotropic magnetization in the vortex 
state of the oxypnictides, particularly the extent to which 
vortex properties are affected by strong magnetic corre- 
lations, multiband effects and possible interband phase 
shift between the order parameters on different pieces of 
the Fermi surface For instance, multiband effects 
in MgB2 can manifest themselves in strong temperature 
and field dependencies for the mass anisotropy parame- 
ter j(T,H) and the London penetration depth X(T,H) 
even at H <C H C 2 ID, 14 1. Yet, there are significant dif- 
ferences between two-band superconductivity in MgB 2 
and in oxypnictides: in MgB 2 the interband coupling is 
weak, while in the oxypnictides it is the strong interband 
coupling which is expected to result in the high T c [B|. 
Thus, probing multiband superconductivity in oxypnic- 
tides by magnetization measurements requires high mag- 
netic fields, which can suppress the superfluid density 
in the band with the largest coherence length above the 



"virtual upper critical field" (H v ) at which the vortex 
cores in this band overlap. In this Letter we address 
these issues, presenting the first high- field torque mea- 
surements of anisotropic reversible magnetization of the 
vortex lattice in SmOo.sFo^FeAs single crystals. Our 
measurements of 7(T, H) up to 30T and extended tem- 
perature range, 20 < T < 40 K have revealed a differ- 
ent behavior of *y(T, H) as compared to recent low- field 
torque measurements [111 ] . 

Measurements of anisotropic equilibrium magnetiza- 
tion m(T, H) in SmOo.sFo^FeAs are complicated by 
the smallness of m(H, T) caused by the large Ginzburg- 
Landau parameter, k = A/£ > 100 and by the strong 
paramagnetism of Sm 3+ ions, which can mask the true 
behavior of m(T,H). In this situation torque magnc- 
tomctry is the most sensitive technique to measure the 
fundamental anisotropy parameters of m(T, H) in small 
single crystals. The torque t = m x H acting upon a 
uniaxial superconductor is given by 
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where V is the sample volume, 0o is the flux quan- 
tum, H$ is the upper critical field along the ab planes, 
rj ~ 1 accounts for the structure of the vortex core, 
is the angle between H and the c-axis, e(0) = (sin 2 9 + 
7 2 cos 2 Q) 1 / 2 and 7 = X c /X a b is the ratio of the London 
penetration depths along the c-axis and the ab-plane. 
The first term in Eq. |T]) was derived by Kogan in the 
London approximation valid at iJ cl <C H <tC U c i [HI- 
The last term in Eq. (1) describes the torque due to para- 
magnetism of the SmO layers and possible intrinsic mag- 
netism of the FeAs layers. Here r m = (x c — Xa)VH 2 /2 
and \c and Xa are the normal state magnetic susceptibil- 
ities of a uniaxial crystal along the c-axis and ab plane, 
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respectively. As will be shown below, the paramagnetic 
term in Eq. (JTJ) in SmOo.sFo^FeAs can be larger than 
the superconducting torque, which makes extraction of 
the equilibrium vortex magnetization rather nontrivial. 
In this Letter we develop a method, which enables us to 
resolve this problem and measure the true angular de- 
pendence of the superconducting torque as a function of 
both field and temperature, probing the concomitant be- 
havior of j(T,H) and X a b{T,H) and manifestations of 
multiband effects in SmOo.sFo^FeAs single crystals. 
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FIG. 1: (color online) (a) Magnetic torque t(8) for a 
SmOo.8Fo.2FeAs single crystal for increasing and decreasing 
angle sweeps (black lines) at 3 T and 27 K. The equilibrium 
T av (8) (blue markers) is obtained by averaging both traces, 
(b) T av (8) for 40 K and 30T exhibits a nearly sinusoidal an- 
gular dependence. Red line corresponds to a fit to the first 
term in Eq. (1) with 7 = 11.5. 



Underdopcd single crystals of SmOi_ x F x FeAs having 
typical sizes, 100 x 100 x 10 ^tm 3 and T c ~ 45K were 
grown by the flux method described in Ref. [16j |. The 
sample was attached to the tip of a piezo-resistive micro- 
cantilever placed into a rotator inserted into a vacuum 
can. The ensemble was placed into a 4 He cryostat cou- 
pled to a resistive 35 T dc magnet of the National High 
Magnetic Field Lab. Changes in the resistance of the 
micro-cantilever associated with its deflection and thus a 
finite magnetic torque r was measured via a Wheatstonc 
resistance bridge. 
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FIG. 2: (color online) (a) Angular dependence of r a v(#) at 3 
T and several temperatures, (b) Same as in (a) for 30 T. 



Fig. 1 (a) shows typical angular dependence of t(6) 
at 27 K and 3 T. Since we are only interested in tem- 
perature and field dependencies of 7 and A, the torque 
data are provided in arbitrary units. A hysteresis, re- 
sulting from the irreversible magnetization is observed 
between increasing and decreasing angle sweeps. Black 
markers depict the average value of both traces, r av (#) = 
(t|(#) + Tj(#))/2 defined as an equilibrium magnetiza- 
tion, where the arrows indicate increasing or decreasing 
angle sweeps. The red line is a fit to the first term in 
Eq. (1) with 7 rs 11.5 a value that is ~ 25% smaller than 



the value reported in Ref. ll|. However, this multipa- 



rameter fit is not very suitable for extraction of the true 
values of 7 due to pronounced error bars for H and a 
significant paramagnetic component particularly at 30T. 
The complete set of the raw T av (8) data is shown in Figs. 
2 (a) and (b). 

The superconducting component of the torque can be 
unambiguously extracted from the data by fitting the 
sum of two measured curves T av (6) + r av (9 + 90°), in 
which the paramagnetic component cancels out: 



t{8) + t(8 + 90°) = 
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FIG. 3: (color online) Angular dependence of T av (#) (blue) 
and T av (# + 90°) (magenta) for 3 T and 30 K. Red line cor- 
responds to a fit of r av (0) + r av (8 + 90°) to Eq. (2). (b) An 
example of the fit of T av (#) to Eq. (1) with the parameters 
taken from Fig. 3 (a). 



where e*(9) = (cos 2 8 + j 2 sin 2 9) 1 / 2 . This procedure is 
illustrated by Fig. 3 (a) where r av (#) for 3T and 30K is 
plotted together with r av (9 + 90°). Black markers depict 
the sum of both traces which is entirely determined by 
the superconducting response. Red line corresponds to a 
fit to Eq. 42]), where r)H^[T] « 315(1 - T 2 /T 2 ) shown 
in the inset of Fig. 4 was extracted from the onset of 
the resistive transition in polycrystalline SmOo.8Fo.2FeAs 
with T c ~ 47 K [l2j|. The onset of the resistive transi- 
tion reflects the behavior of those crystallites having the 
field along the ab-plane thus i/° 2 b (T). Given the lack of 
high-field H c2 (T) data for our single crystals, the use of 
the measured value for polycrystals eliminates the 
ambiguities of the multiparameter fit if H C 2 in Eq. (2) 
is treated as one more fit parameter. Using the super- 
conducting parameters obtained by this method in Eq. 
(1) we then fit the original T av (9) (red line in Fig. 3 
(b)), treating r m as the only adjustable parameter. This 
method does result in an excellent fit, allowing us to ex- 
tract both the superconducting torque and the param- 
agnetic torque associated with Sm or Fe moments. It 
also gives smaller values for 7(T) than the ones obtained 
by direct fitting r av (9) uniquely to the first term of Eq. 
(fT]). As shown in Figs. 4 and 5 (a), the so-obtained 7 is 
not only temperature dependent but it is also field 
dependent, decreasing by more than 20 % at 30 T. 

The obtained temperature dependence of j(T,H) is 
reminiscent of the behavior of j(T, H) previously re- 
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FIG. 4: (color online) Temperature dependencies of 7 ob- 
tained from fitting T a v(d)+T av (8+90°) to Eq. (2) for 3T (solid 
circles) and 30T (open circles). The ambiguity of the extrac- 
tion of 7 from a direct multiparameter fit of r av (6) [ll|] to the 
first term in Eq. (1) is illustrated by the triangles. Inset shows 
the difference of H C 2 for different fit procedures. Squares cor- 
respond to r]Hc 2 obtained by fitting the raw r a v(8, T, H — 3T) 
solely to first term in Eq. (I). Solid circles correspond to H^2 
measured on a polycrystalline SmOo.8Fo.2FeAs with a T c ~ 47 
K [TJ| which we used to fit the data to Eq. (2) and from which 
we obtain H^ 2 = /j (open circles are points from H = 3 
T while crossed circles from H = 30 T). 



ported for MgB2 which was explained in terms of multi- 
band effects [131,1141 • Yet the extracted London penetra- 



tion depth shown in Fig. 5 (a) does not exhibit a signif- 
icant field dependence which would indicate an abrupt 
depression of the superfluid density in one of the bands 
above H v (produced by the suppression of the respective 
superconducting gap). This is quite remarkable given 
that we measured A(T, H) up to the applicability limit of 
the London theory, i.e. r/H^/H ~ 1 at H = 30 T. Over- 
all, the behaviors of 7(T, H) and A(T, H) shown in Fig. 
5 (a) would suggest two strongly coupled gaps of similar 
magnitude but not too different mass anisotropics. The 
decrease of 7(H) as H increases may indicate that the 
band with the shorter coherence length £ ~ Hvf/A is 
the least anisotropic. 

The relative contributions of the superconducting and 
magnetic components in T av (0, T, H) are shown in Fig. 5 
(b). At higher T the behavior of r m (T) oc d/T[K] + C 2 , 
C2 ~ — Ci/43 at 3T is consistent with the Curie- Weiss 
paramagnetism of Sm 3+ ions. However, this tempera- 
ture dependence of r m (T) changes at H = 30 T, for 
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FIG. 5: (color online) (a) Field dependence of 7 (blue mark- 
ers) and C\ a b (green markers) obtained by fitting (t(0)+t(6+ 
90°) at T = 27 K to Eq. (2). Here, C is a calibration de- 
pendent constant for our cantilever. Black markers depict the 
field dependence of respectively the magnetic (open squares) 
and superconducting (solid squares) components. For clar- 
ity the data was multiplied by a factor of two. Red line 
corresponds to a T 2 fit and the magenta line to a T 3 de- 
pendence, (b) Amplitudes of superconducting (circles) and 
magnetic (triangles) components of r av (6), at H = 3 T (blue 
markers) and 30 T (red markers). 



where n is the density of paramagnetic ions, gg = 
(gl cos 2 9 + gl sin 2 9) 1/2 . For weak fields or fisgeH < T, 
Eq. © gives t p = (xc - Xa)H 2 sin 20/2 and Xa = 
figg^/AT as used in our analysis. However, for higher 
fields fisgeH > T, the paramagnetic torque t p ~ 
n^s(g 2 — g%)H sin 26/2gg acquires higher order harmon- 
ics. This case may pertain to our data at 30T and 
T < 30K, for which the pure sin 29 component in t(9) 
does decrease as T decreases, but t p (9) may not be com- 
pletely eliminated by the procedure described above. De- 
viations from Eq. |T]) also come from corrections to the 
London theory at high fields resulting in additional terms 
oc aH c2 /H - (Inr) + a)H/H c2 , a ~ 1 in m(T, H, 9) due 
to pairbreaking and nonlocal effects [2(| . 

In summary, our torque measurements at high-fields 
reveal the temperature and field dependencies of the 
anisotropic reversible magnetization which is strongly 
coupled with the magnetism of rare earth ions in 
SmOo.sFo^FeAs single crystals. Our results indicate 
a temperature and field dependent mass anisotropy 
7(T, H) which saturates at 7 ~ 9 at low temperatures 
under a modest field. This value is higher than 7 = 
H^ 2 /H^2 — 5 at low temperatures in NdO .7F . 3 FeAs 
single crystals [2l| and 7 ~ 5 — 7 for YBa 2 Cu3Q7 $, but 
is much smaller than 7 > 30 suggested by Ref. [22j . The 
observed insensitivity of the London penetration depth at 
fields up to 30T is indicative of strong coupling supercon- 
ductivity, which in addition to a not very high 7 is very 
important for applications. Our results are consistent 
with strongly coupled gaps of nearly equal magnitudes 
in distinct bands. 
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which T m (T) decreases as T decreases, indicating the in- 
triguing possibility of field-induced antifcrromagnetism 
(AF). Indeed, the coexistence of AF and superconduc- 
tivity in phase-separated regions has been reported in 
the oxypnictides [171 ] . while AF in the vortex cores of 
the cuprates has been observed at higher field L8(. In 
our case, the latter can be ruled out since the ob- 
served T m (H,T = 27K) follows a H 2 dependence, un- 
like r m oc H 3 /H C 2 for the AF cores. Another mecha- 
nism may result from a field-dependent uniaxial magnetic 
anisotropy as the Sm 3+ moments align along the field. 
This effect can be modeled by the single-ion Hamilto- 
nian H = HB(g a cr x H x + g a a v H y + g c a z H z )/2, where a a 
are the Pauli matrices, and g n are the principal values of 
the effective g-factor tensor [19| . In this case 

t p = nfi B {g 2 c - g 2 a )Hsm29t&nh{ti B g e H/2T)/2gg, (3) 
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